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1 Introduction 

When the size of materials is comparable to the characteristic length scale of their physical 
properties, novel functionalities can emerge. For semiconductors, this is exemplified by 
the "superlattice" concept of Esaki and Tsu, where the width of the repeated stacking of 
different semiconductors is comparable to the "size" of the electrons, resulting in novel 



confined states now routinely used in opto-electronics 



jj. For metals, a good example is 



magnetic/non- magnetic multilayer films that are thinner than the spin-scattering length, 
from which giant magnetoresistance (GMR) emerged 2|, |3|, used in the read heads of 
hard disk drives. For transition metal oxides, a similar research program is currently 
underway, broadly motivated by the vast array of physical properties that they host. This 
long-standing notion has been recently invigorated by the development of atomic-scale 
growth and probe techniques, which enables the study of complex oxide heterostructures 
approaching the precision idealized in Fig. dl^a). Taking the subset of oxides derived from 
the perovskite crystal structure, the close lattice match across many transition metal 
oxides presents the opportunity, in principle, to develop a "universal" heteroepitaxial 
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Figure 1: (a) Schematic illustration of ideal heterointerfaces of two perovskites ABO3 and 
A'B'Oa stacked in the [001] direction, (b) Charge sequences of the AO and BO2 planes 
of perovskites plotted together with their pseudocubic lattice parameters, (c) Scanning 
transmission electron microscopy image of a LaTiOa/SrTiOa (001) superlattice (Ohtomo 
et al. fl). 



materials system. 

Hand-in-hand with the continual improvements in materials control, an increasingly 
relevant challenge is to understand the consequences of the electrostatic boundary con- 
ditions which arise in these structures. The essence of this issue can be seen in Fig. HJ^b), 
where the charge sequence of the sublayer "stacks" for various representative perovskites 
is shown in the ionic limit, in the (001) direction. To truly "universally" incorporate 
different properties using different materials components, be it magnetism, ferroelectric- 
ity, superconductivity, etc., it is necessary to access and join different charge sequences, 
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labelled here in analogy to the designations "group IV, III-V, II- VI" for semiconductors. 
As we will review, interfaces between different families creates a host of electrostatic 
issues. They can be somewhat avoided if, as in many semiconductor heterostructures, 
only one family is used, with small perturbations (such as n-type or p-type doping) 
around them^. However, for most transition metal oxides, this is greatly restrictive. For 
example, LaMnOa and SrMnOs are both insulators in part due to strong electron corre- 
ations, and only in their solid solution does "colossal magnetoresistance" emerge in bulk 
6| . Similarly, the metallic superlattice shown in Fig. [T]^c) can be considered a nanoscale 
deconstruction of (La,Sr)Ti03 to the insulating parent compounds. Therefore the aspira- 
tion to arbitrarily mix and match perovskite components requires a basic understanding 
of, and ultimately control over, these issues. 

In this context, here we present basic electrostatic features that arise in oxide het- 
erostructures which vary the ionic charge stacking sequence. In close relation to the 
analysis of the stability of polar surfaces and semiconductor heterointerfaces, the vari- 
ation of the dipole moment across a heterointerface plays a key role in determining its 
stability. Different self-consistent assignments of the unit cell are presented, allowing 
the polar discontinuity picture to be recast in terms of an equivalent local charge neu- 
trality picture. The latter is helpful in providing a common framework with which to 
discuss electronic reconstructions, local-bonding considerations, crystalline defects, and 
lattice polarization on an equal footing, all of which are the subject of extensive current 
investigation. 



2 The polar discontinuity picture 
2.1 Stability of ionic crystal surfaces 

The surface of crystals determines many of their physical, mechanical and chemical prop- 
erties. Due to the lack of translational symmetry in the perpendicular direction, the 

^These effects can in principle also be reduced by choosing a (110) growth orientation, but other 
aspects of stability may be limiting Q . 
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Figure 2: Schematic illustration of the charge density p, the dipole fi in the unit cell 
starting from the top-most layer, the electric field E induced by the dipoles, and the 
electrostatic potential (p at the surface of an ionic crystal with dipole moment in each 
unit cell. 



stable charge distribution at the surface can be completely different from that of the 
bulk, and the surface may reconstruct in a manner different from the bulk states. Imag- 
ine an ideal ionic crystal which consists of charged ions bound together by their attractive 
interactions, and all the ions are taken as fixed point charges. Since the charges are lo- 
cally preassigned to the ions in this model, the ideal ionic surface apparently requires no 
reassignment of the charges from that of the bulk. However, the electrostatic potential 
in an ionic crystal diverges when there is a dipole moment in the unit cell perpendicular 
to the surface. The potential^ (p should be constant in vacuum in the absence of external 
fields, and the potential can be obtained by integrating the electric field caused by the 
charged sheets, as shown in Fig. [2l A finite shift in the potential emerges due to the 
dipole moment of each unit cell, and as the unit cells are stacked, so the potential grows, 
and diverges into the crystal. Due to this effectively infinite surface energy, such surfaces 
cannot exist without reconstructions, and the stability of an ionic surface randomly cut 
from the bulk is not trivial without knowing the stacking sequence of the charged sheets 
precisely. This surface instability and the associated reconstructions have been indeed 
observed by means of low-energy electron diffraction (LEED) and ion scattering, where 
absorption of foreign atoms, surface roughening, or changes in surface stoichiometry were 
found \m\- 

In order to survey the stability of such surfaces, Tasker introduced a classification of 



^Notc that literature on this topic uses both the electrostatic potential (for a positive test charge) 
and the electron potential energy (as in band diagrams) — we use the former here. 



the surfaces from the viewpoint of the charge o: 
unit cell which starts from the top-most layer 



the atomic sheets and the dipole in a 



10|. Note he discussed the stability of bulk 



frozen surfaces'^, where the top-most layer is one of the constituent atomic sheets of the 
bulk crystal and has no reconstruction. Tasker described three types of the surfaces, as 
shown in Fig. [3l 

• Type 1 has equal numbers of anions and cations on each plane, and therefore the 
unit cell has no dipole moment. For example, the (001) and (110) surfaces of the 
rocksalt structure MX (e.g. NaCl, MgO, NiO) are classified as this type. 

• Type 2 has charged planes, but no net dipole moment perpendicular to the surface. 
The anion X terminated (001) surface of the fiuorite structure MX2 (e.g. UO2, 
Th02) is an example. 

• Type 3 has charged planes and a net dipole moment normal to the surface. Exam- 
ples include the (111) surface of the rocksalt structure, and (001) or (111) surfaces 
of the zincblende structure MX (e.g. GaAs, ZnS). 

In Tasker 's classification, type 1 and type 2 surfaces are stable while type 3 is not, 
since the instability of the surface comes from the stacking of the dipole in each unit 
cell. Macroscopically, the instability of the type 3 surface arises from the change in 
the potential slope when crossing the surface. The term "polar surface" can be defined 
following this classification, namely we call a surface "non-polar" when it is type 1 or 
type 2, and "polar" when it is type 3. 

These definitions are analogous to the definition of polar crystals, although polar sur- 
faces and surfaces of polar crystals are not equivalent. Dielectric polarization is observed 
when an electric field is applied to a material, but even in the absence of the field, some 
crystals retain a "spontaneous" polarization 12|. Only 10 out of 32 point groups show 



this behavior, and their members are called polar crystals, while the others are non-polar. 
Twenty-one out of 32 point groups do not have inversion centers, and the polar crystals 



This definition of ^'bulk frozen" follows description of Goniakowski et al. [11 
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Figure 3: Distribution of charges p on planes for three stacking sequences parallel to the 
surface, (a) Type 1, (b) type 2, and (c) type 3 (Tasker (lol|). 



are included among them. When a material shows macroscopic spontaneous polariza- 
tion, the electrostatic potential of one end is different from the other end as a result of 
the stacking of the dipole in each unit cell, and the surface usually has "compensating" 
charge to reconcile this potential difference. Since Tasker took unit cells from the top- 
most layer, even crystals with inversion symmetry can show dipoles in the unit cells in 
his model. For example, although NaCl is cubic and has inversion symmetry, its (111) 
surface is classified into type 3. Therefore, the word "polar" should be used with some 
care since it has different meanings in different contexts. 
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2.2 Stability of covalent surfaces 

At the surface of covalent crystals, lacking full coordination, the top-most atoms have 
valence electrons which are not used for bond formation. These non-bonding electrons 
are called dangling bonds, and have higher energy than the bonding electrons, which 
causes the movement of the atom positions to decrease their number 

In a covalent crystal, since the bonds are formed as a hybridization of the valence 
electrons of charge-neutral atoms, one can describe the charge distribution starting from a 
bulk unit cell which is charge neutral and dipole-free. However, when the electronegativity 
of the atoms are different e.g., Ga and As in GaAs, charge transfer between anions and 
cations occurs, similar to the ionic case. This charge transfer is realized by the difference 
of the contribution of each bond, namely 1 + a electrons to the anions and 1 — a electrons 
to the cations. Here a is a parameter to describe the ionicity of the bond, determined by 



the electronegativity of the two bonding atoms. As a consequence, the unit cell can have 
a dipole moment normal to the surface, which causes the same instability as that in the 
ionic crystal case. Therefore, even in a covalent crystal, a surface instability emerges from 
the dipole in the unit cell, independent of the surface instability naturally arising from 
dangling bond formation. Even though Tasker's classification was introduced to describe 
the stability of ionic surfaces, it is also relevant for covalent surfaces in the presence of 
finite ionicity. 

Both dangling bond formation and the instability of polar surfaces are at play, and 
they are reconciled simultaneously at the surface of covalent crystals. Therefore, direct 
observation of the instability of polar surfaces in covalent systems has been difficult. 
When we consider an epitaxial interface which has similar charge structure as the polar 
surface, the instability from the dangling bonds disappears, and we can solely discuss the 
stability in the same manner as that for the ideal ionic surfaces, as discussed in the next 
section. 



2.3 Polar semiconductor interfaces 

Similar to the instability of polar surfaces, dipoles in the unit cells stacking from the 
interface can cause potential divergence and instability, and require reconstruction. This 
point was first proposed at the heteroepitaxial interface of GaAs and Ge in the [001] 
direction by three groups from different starting points for treating ionicity. Based on 
their considerations, we define polar and non-polar interfaces, in analogy to polar and 
non-polar surfaces. 



2.3.1 Charge transfer based on electronegativity 

Frensley and Kroemer calculated the energy band diagram at abrupt semiconductor het- 
erojunctions ij]. Their starting point was to describe the alignment of atoms around the 
interface without considering the ionicity, and then calculate the charge transfer based on 
the ionicity using the electronegativity of the atoms, under assumption that the charge 
transfer only occurs between the nearest neighbors. The Phillips electronegativity values 




Figure 4: (a) Model for a (001) Ge/GaAs heterojunction considering the ionic character 
of the bonds. The atomic positions and effective ionic charges are shown above. Below 
is a diagram of the plane-averaged potential (Frensley and Kroemer (b) Schematic 

diagram of the charge transfer from the neutral atoms with respect to the electronega- 
tivity. 



Xph were used (Xph(Ga) = 1.13, Xph(Ge) = 1.35, and Xph(As) = 1.57 (isf). In the bulk 
zincblende structure AB, the A site is tetrahedrally coordinated by four B atoms and 
vice versa, and based on their calculation 



16| , the ionic charges e* of the atoms are given 



by 



e*iA) = -e\B) = 0.76go x [X^^{B) - Xph(A)], 



where go is the elementary charge. This is equivalent to assuming a charge transfer of 
0.76go X ^[Xph(S) - Xph(A)] between any pair of the nearest neighbors. 

Consider the case of Ge/GaAs interfaces as shown in Fig. Ufa), where the charge 
e*(Gaint) on the Ga ions adjacent to the interface is 



e*(Gaint) = 0.76go x 



^Xph(Ge) + ^Xph(As) - Xph(Ga) 



0.25go. 



(2) 



Similarly, the charges e*(Geint) on the Ge ions at the interface and e*(As) at the As sites 
are e*(Geint) = -O.OSgo and e*(As) = -0.33go- 

Without ionicity, the electrostatic potential is constant, and even with ionicity, the 
potential does not diverge. This can be understood easily by tracking the virtual charge 
transfer processes from the starting alignment of the charge neutral atoms. The charge 
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Figure 5: (a) Calculated contour plot of charge density for the partially occupied interface 
band around Ga-Ge bonds (Baraff et al. (b) Schematic model for counting electrons 

from the local-bond point of view. 

transfer is equivalent to the situation that each neutral atom loses 0.76go x -^Ph charges, 
and half of them are transferred to the left atoms and the other half to the other side, 
as shown in Fig. ID^b). Therefore, charges are always transferred symmetrically, and each 
modulation creates no dipole, resulting in no potential shift. Here, a change of the number 
of electrons in the ions compared to the bulk state is assumed, which is equivalent to 
changing the valence assignments. 

2.3.2 Self-consistent calculation and counting electrons of bonds 

Baraff et al. performed a self-consistent calculation of the potential, charge density, and 
interface states for the abrupt interface between Ge and GaAs, terminated on a (001) 

n 

Ga plane [17|. As shown in Fig. El^a), their calculation showed fractional occupancy of 
electronic states only at the interface, which cannot exist in bulk. 

These interface states can be discussed from a local-bond counting point of view as 
well. The number of the valence electrons is 4 for Ge, 3 for Ga, and 5 for As. When we 
assume all the valence electrons of an atom are equally distributed to the four covalent 
bonds around it, Ge, Ga, and As atoms supply 1, 0.75, and 1.25 electrons to each bond, 
respectively. As shown in Fig. |5]J^b), the Ge-Ga bonds at the interface have 1.75 electrons 
per bond, while other bonds have 2 electrons in each. These partially occupied bonds 
are considered to form the interface states. 

When the number of electrons are smaller than that in the bulk, the attractive inter- 
action between the bonded atoms should be weaker. According to their calculation, the 
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energy is minimized when the Ge-Ga bond length is 4 % larger than that of the bulk. 
Without the change of the bonding length, the system requires long-range disturbance 
of the lattice, which is unlikely to be realized. 



2.3.3 Solving the potential divergence from the ionic picture 

Both the charge transfer model by Frensley and Kroemer and the electron counting model 
by Baraff et al. predict (indeed require) that there are interface states at a Ge/GaAs (001) 
interface even if it is perfectly abrupt, with no crystalline defects. The central point raised 
by these studies is that despite having the same crystal structure, and having very closely 
matching lattice constants, this interface must accomodate charge arising from interface 
boundary conditions. However, experimentally no considerable density of interface states 



was observed 18|, and a model to treat this interface without changing the number of 
charges at the interface was required. 

Harrison et al. constructed a model for the Ge/GaAs (001) heterojunctions by stack- 
ing the fully ionized atoms, and calculated the electrostatic potential based on the fixed 
assignment of the charges 



19|. As shown in Fig. |6]^a), the potential is very similar to the 
case of polar surfaces since the unit cells which start from the interface have dipoles in 
GaAs, while the unit cells in Ge are always charge neutral. Therefore, the stacking of 
the dipoles causes potential divergence in this case as well. 

Since the charge of each ion was fixed, the solution to the instability of this interface 
requires compensation by changing the stoichiometry at the interface. They proposed a 
simple model where 1/4 of the Ge atoms are replaced by As atoms at the interface while 
1/4 of the Ga atoms adjacent to the interface are replaced by Ge. In this reconstructed 
model with two transition layers, the electrostatic potential does not diverge, as shown 
in Fig. El^b). 

It might be surprising that from two completely different starting points, namely one 
from covalent (Frensley and Kroemer) and the other from ionic (Harrison et al.) pictures, 
exactly the same potential diagrams were obtained. However, rearranging the number 
of charges at the perfectly abrupt interface or changing the interface composition while 
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Figure 6: Schematic crystal structure and electrostatic potential (f) in the heteroj unctions 
of Ge and GaAs in the [001] direction, (a) An atomically abrupt interface, (b) A Ge/GaAs 
heteroj unction with two off-stoichiometric transition layers (Harrison et al. [l9|). 

maintaining the ionic charges of the atoms can give the same net charge distribution. The 
experimental absence of localized states suggests the atomic reconstruction based on the 
ionic picture. This is equivalent of saying that the electronic state at this semiconductor 
heterointerface cannot deviate so strongly from that of the bulk constituents — it is en- 
ergetically inaccessible. This is the fundamental aspect which can be quite different in 
complex oxide heterointerfaces, and is the subject of much current excitement. Namely, 
there is a possibility that the charge transfer picture (Frensley and Kroemer) and the 
electron counting picture (Baraff et al), which require large deviations of electron num- 
bers from the bulk values, can be energetically accessed in transition metal oxides with 
multi- valency, as described in Section |31 

2.3.4 Definition of the polar interface 

Following the model by Harrison et al, we can define the polar nature of a bulk frozen 

interface between two materials, where the interface consists of two bulk frozen surfaces 

connected together^. First let us classify the interfaces by the polar nature of the two 

constituent surfaces, as shown in Fig. [71 

*Here, we consider only interfaces between two scmi-infinitc materials — we ignore the coupling to 
other interfaces or surfaces, which is discussed in Sections 14.31 and 14.41 
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Figure 7: Schematic distribution of charges p on planes, the dipole moment fj, in unit 
cells starting from the interface, and the electrostatic potential for the four stacking 
sequences parallel to the interface, (a) type I, (b) type II, (c) type III, and (d) type IV. 
In order to treat a bulk frozen interface as a set of two bulk frozen surfaces, vacuum is 
inserted between them (dashed lines). The electrostatic potential was calculated taking 
the vacuum as the potential reference except for the right plot in (c), where a constant 
electric field was added to show the absence of macroscopic band bending at the interface. 



• Type I is formed between two non-polar surfaces, and the potential is flat. 

• Type II is formed between polar and non-polar surfaces, and the potential diverges 
from the interface in the material with the polar surface. 

• Type III is formed between two polar surfaces, and the direction and the size of the 
dipoles in the unit cells which start from the interface is the same. Although (p looks 
to diverge from the interface in both materials, there is no macroscopic difference 
in the potential slope across the interface. We can cancel the potential slope on 
both sides by adding a constant electric field, as shown in Fig. El^c). Therefore, this 
interface is stable as constructed. 

• Type IV is formed between two polar surfaces, where the dipoles in the two different 
unit cells are not identical, which results in a macroscopic difference in the potential 
slope at the interface. It is impossible to find any constant electric field to cancel 
out the potential divergence in both media, due to this difference^. 

■'Note type IV is the most common and general case in reality, since the electronegativity can never 
perfectly match between different materials. 
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Type I and type III are stable due to the absence of a macroscopic difference of the 
potential slopes, while type II and type IV are not stable. This difference arises from the 
continuity/ discontinuity of the dipoles in the unit cells which start from the interface. 
A bulk frozen interface is non-polar, when the dipoles in the unit cells starting from 
the interface are identical across the interface, and thus no change of the macroscopic 
potential slope exists. On the other hand, it is polar if a bulk frozen interface has a 
discontinuity of the dipole moment in each unit cell. For example, from a purely ionic 
viewpoint (Harrison et al.) the abrupt Ge/GaAs (001) interface is classified as type II, 
and thus polar. 

This definition of the polar nature of interfaces is consistent with that of surfaces. 
When the vacuum is treated as a charge neutral medium, the non-polar surface is type 
I, and the polar surface is type II in this classification of the interfaces, and the polar 
nature is maintained following the definitions for the interface. We can treat surfaces 
and interfaces in one framework, which is the polar nature of discontinuities at materi- 
als boundaries. In summary, polar discontinuities, which consist of polar surfaces and 
interfaces, are unstable due to the macroscopic potential folding arising from the discon- 
tinuity of the stacking of dipoles in the unit cells, and require reconstructions to stabilize 
them. Interfaces with continuous polarity, on the other hand, are stable without any 
reconstructions. 



3 Metallic conductivity between two insulators 

As noted in the introduction, the modern ability to approach atomic control in complex 
oxide heterostructures has newly enabled the experimental investigation of their polar 
discontinuities. As illustrated in Fig. [D^a), an immediate question arises regarding the 
choice of the termination layer at the interface, and the consequences of this degree of free- 
dom. This issue has been most explicitly raised, and hotly debated, for the electron gas 



observed at the interface of two perovskite insulators, LaAlOs and SrTiOs 20| . Specif- 
ically, the (001) heterointerface was found to be insulating when grown using a SrTiOs 
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substrate which was SrO-terminated, and conducting when Ti02-terminated. Given the 
rapid evolution of the field, and the numerous reviews of this heterostructure in the lit- 
erature, we do not attempt a comprehensive review here. Instead the LaAlOa/SrTiOs 
interface will be used to illustrate the various mechanisms suggested to explain the in- 
terface electronic structure, and the electrostatic boundary conditions which arise. It 
should be stressed, however, that all oxide heterointerfaces should be considered type IV 
to varying degrees, and thus these issues are quite general - even underlying the interface 
charge in the superlattice shown in Fig. [U^c). 



3.1 The polar discontinuity scenario 

Assuming pure ionicity, the charge sequence of the (001) perovskite planes are different 
in these two materials, namely the planes of LaAlOs are (La^"''0^~)"^ and {Al'^^Ol~)~ , 
while those of SrTiOa are (Sr^+O^")^ and (Ti^+Og")'^. Therefore, the abrupt interface 
between LaAlOa and SrTiOa is type II polar^ and requires reconstruction as shown in 
Fig. [HI just as for the GaAs/Ge (001) interface. 

Unlike polar semiconductor interfaces where only atomic reconstructions are avail- 
able due to the fixed ionic charge of each element'', we have another possibility to rec- 
oncile the instability of polar interfaces, through electronic reconstructions 23|. At the 
LaAlOs/SrTiOs interface, when the interface termination is LaO/Ti02, it requires a net 
half negative charge per 2D unit cell to reconcile the potential divergence (n-type). Ac- 
cessing Ti^"*" can source this negative charge by accommodating electrons at the Ti 3d 



level, as was spectroscopically observed 



24l |. On the other hand, when the interface is 



terminated by A102/SrO, the sign of the required charges is opposite (p-type). Due to 
the difficulty of accommodating holes in this structure (such as Ti^"*"), the positive charges 
are realized by the formation of oxygen vacancies close to the interface, as inferred from 
measurements of the 0-K edge fine structure. Oxygen vacancies are known as electron 



^Allowing for covalency, SrO and Ti02 planes in SrTiOs arc no longer charge neutral, and thus the 
(001) surface is weakly polar, but still the dipolc size of the unit cells which start from the interface is 
different from that of LaAlOs. 

^While here we discuss large scale charge modifications, small polar discontinuities can induce free 
carriers in semiconductors, such as in AlGaN/GaN heterostructures pll [2^. 
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Figure 8: Polar reconstructions at the LaAlOa/SrTiOs interfaces. The unreconstructed 
(a) LaO/Ti02 terminated interface, (b) A102-SrO terminated interface, (c) and (d) the 
corresponding reconstructed interfaces, respectively (Nakagawa et al. |2J|). 



donors, but in this case they are formed to provide positive charges, and thus no electrons 
are supplied from these vacancies compensating the polar discontin uity . Therefore the 
system does not have itinerant electrons and remains insulating 



20 
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One of the key corollaries of this scenario is the LaAlOs thickness dependence. This 
is because the size of the potential difference arising from the stacking of the dipoles in 
the unreconstructed structure is finite in thin films, and if it is small enough, the system 
may be stable without any reconstruction. Indeed, a critical thickness of LaAlOa was 
observed 



261], where the n-type LaAlOs/SrTiOa interface is insulating if the thickness of 
LaAlOs is up to 3 unit cells, and metallic above that thickness. A similar tendency was 
also observed in SrTiOa/LaAlOa/SrTiOs trilayer structures where the distance between 



the two polar interfaces was varied 



27|. 



3.2 Oxygen vacancy formation during growth 

SrTiOs is known to be a material which readily accommodates oxygen vacancies that 



act as donors to provide itinerant electrons 
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29|. Either kinetic bombardment of the 



SrTiOs substrate by the ablated species (early studies of this interface all used pulsed 
laser deposition), or gettering by a reduced film, can induce oxygen vacancies, and they 
were suggested to be the dominant origin for the observed conductivity by several groups 
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Figure 9: Conducting tip atomic force microscopy mapping around the LaAlOs/SrTiOs 
interface of (a) the as grown sample {P02 = 10~^ Torr) and (b) the postannealed sample 
(Basletic et al. H). 



30l-l32|. Indeed, the first report found a strong dependence of the Hall density for n-type 
interfaces on the oxygen partial pressure (-P02) during growth, while the p-type interface 



was robustly insulating (20|. For n-type samples with similar variations in the transport 
properties, Kalabukhov et al. found when grown at P02 = 10^^ Torr, the samples 



exhibited blue cathode- and photo-luminescence at room temperature 



30|, similar to 



that of reduced SrTiOa by Ar bombardment 33j. In addition to transport studies, 
Siemons et al. demonstrated that the photoemission spectra from these interfaces showed 
a larger amount of Ti^"*" in samples grown at low pressures without oxygen annealing 



3l| . Herranz et al. observed Shubnikov-de Haas oscillations in reduced LaAlOa/SrTiOs 
samples which were quite similar to bulk doped SrTiOa, and rotation studies indicated a 
three-dimensional Fermi surface 32 1. 



The strong P02 dependence of the conducting channels in LaAlOs/SrTiOs was ob- 
served and mapped by means of conducting tip atomic force microscopy on cross-sections 
of the interface, which revealed a conducting region extending > 1 //m into the substrate 



for samples grown at low pressure 



34l |. as shown in Fig. [HI After annealing, the width 



of the conductive layer decreased to ~ 7 nm, as limited by the radius of the probing tip. 
The consensus of these and other studies was that the free carriers in samples grown at 
low P02 were dominated by oxygen vacancies, since the density far exceeded that needed 
to stabilize the polar discontinuity. For high P02) after post-annealing, the origin was 
less clear. 
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Figure 10: (a) Occupancies and (b) cumulative displacements of the atoms at the 
LaAlOs/SrTiOs interface, (c) Concentration of Ti^+ determined by a minimization of 
the electrostatic potential, (d) Predicted cumulative unit cell displacements from bulk 
positions, based on the component ionic radii (Willmott et al. Q). 



3.3 Intermixing and local bonding at the interface 

Willmott et al. studied a five unit cell film of LaAlOa on SrTiOa (00 
surface x-ray scattering 



3y means of 



36| and further 



351 ] ■ using coherent Bragg rod analysis (COBRA) 
structural refinement. Their analysis revealed intermixing of the cations (La, Sr, Al, 
and Ti) at the interface [Fig. nUT a)]. as well as significant local displacement of the 
atomic position both in the film, and in SrTiOs close to the interface [Fig. [TOT b)]. The 
distribution of Ti valence was also inferred by minimizing the electrostatic potential 
[Fig. [TUT c)] following the obtained atomic positions, and the atomic displacements were 
explained by the larger ionic radii of Ti'^"'" compared to Ti^"*" [Fig. [TOT d)]. Based on these 
observations, the origin of the interface conductivity was suggested to be the formation of 
the bulk-like solid solution Lai-ajSr^-TiOs in a region of approximately 3 unit cells. This 
explanation can be considered a diffused version of local bonding arguments — i.e., that 
even in the abrupt limit, the Ti at the interface has La on one side, and Sr on the other. 



3.4 Reconciling the various mechanisms 

At present, it can be fairly stated (we believe) that no one scenario can completely ex- 
plain the vast and growing body of experimental work on this system. Even theoretical 
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calculations show an extreme sensitivity to the choice of boundary conditions and as- 
sumptions of site-occupancy 
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4l| . To give examples for each perspective: The polar 



discontinuity 
experiments 



picture should lead to a significant internal field in ultrathin LaAlOs, while 



42| put an upper bound far below that expected theoretically, even allowing 



for significant lattice polarization 



43| . Oxygen vacancies induced by growth are difficult 



to reconcile with the notion that a single layer of SrO can prevent their formation. Local 
bonding and interdiffusion considerations do not address the constraints of global charge 
neutrality. Furthermore, discriminating between these mechanisms is often difficult, since 
the change of Ti valence shows similar transport, spectroscopic, and optical properties, 
independent of its origin. 

It is extremely likely that multiple contributions exist to varying degrees, dependent 
on the growth details of a given sample. While this is a matter for further experimental 
investigation and refinement, an equally difficult issue appears to be one of semantics. 



For example, one of the conceptual difficulties of the polar discontinuity picture 
been the question: Where do the electrons come from? In many presentations 
26, 



has 
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45| . the charges at the interface are depicted to originate from the surface of 



the LaAlOs film, but it is not so obvious that they travel through the insulating films 



independent of its thickness 



261]. Fundamentally, the "non- locality" of these electrostatic 



descriptions has sometimes been considered less intuitive and compelling than "local 



chemistry" mechanisms such as vacancies or interdiffusion 



30, 



31 



34 



35, 



46| 



A related difficulty is how to treat the charge density to describe the macroscopic 
electric field. Based on the polar discontinuity picture, the stacking of the dipoles in 
the unit cells creates a macroscopic electric field, resulting in the change of the potential 
slope at the discontinuity. In this picture, the unit cells start from the discontinuity, and 
thus the composition is always the same as that of the bulk, which is charge neutral. It 
might be strange that the potential starts to bend at the discontinuity although all the 
unit cells start out charge neutral, because the source of an electric field is charge. In 
fact, the source is the macroscopic bound charge density at the interface, as discussed 
in Section but the existence of such implicit charge density has caused a fair bit of 
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confusion. To address these concerns, it is useful to treat the electrostatics in a purely 
local description, as well as the boundary charges explicitly. Furthermore, the effects of 
defects and diffusion can be discussed more simply by re-framing the polar discontinuity 
picture in local form. Therefore, this local charge neutrality picture based on dipole-free 
unit cells is developed in the next section first for idealized models, followed by discussion 
of incorporation of chemical defects. 



4 The local charge neutrality picture 
4.1 Unit cells in ionic crystals 

One of the origins of confusion regarding the stability and reconstructions of polar discon- 
tinuities arises because of the various choices of a unit cell in a crystal, which determines 
the size and direction of the dipole in it^. For example, when we take the unit cell of 
LaAlOa, this stacking can be treated as a dipole of [(A102)^ (LaO)''"], as shown in Fig. 
[nja). However, it is also possible to take [(LaO)"^ (A102)^] as a unit cell, and the sign 
of the dipole is opposite to the previous case, as shown in FigJTTlfb). 

The choice of unit cells should not change the electrostatic potential in the crystal. 
Indeed, the difference between the two choices is compensated by the potential 0sur arising 
from the surface layer. If the surface is terminated by the (A102)~ layer, it remains as an 
extra negatively charged layer when we take [(LaO)"^ (A102)^] as the unit cell, as shown 
in FigJTlTb). and the total electrostatic potential remains the same as in the [(A102)~ 
(LaO)"'"] unit cell case. 

Therefore, it is impossible to fix the direction of the dipoles without knowing the 
surface termination. In other words, it is the surface that determines the stability of 
the system. So the problem can be simplified by considering the surface locally, and not 
by counting the number of dipoles in the material. In order to avoid the dipoles arising 



^For simplicity, we started our diseussion from a point-eharge model of ions, and neglect the free 
carrier distribution. More generally, the charge distribution can be described using contributions from 
ion cores and free carrier Wannier functions, and the ambiguity of the choice of unit cells appears in this 



extended case as well [47 1. 
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Figure 11: Schematic illustrations of the various choices of unit cells in LaAlOs, charge 
distribution p, and associated potential 0. The filled allows indicate the orientation of 
the dipoles with size A in the unit cells, (a) Taking [(A102)~ (LaO)"'"] as a unit cell, 
(b) Taking [(LaO)"'" (A102)~] as a unit cell, and the total potential (ptot is the sum of 
the potential (pdip arising from the stacking of the dipoles and the potential (psm from the 
surface charge, (c) Taking a dipole-free unit cell [|(LaO) - (AIO2) - |(LaO)]. 




le simplest approach is to take dipole-free unit cells. 



from the stacking of charged layers, t 
This can be achieved in any crystal and in the LaAlOs case, this is done by taking 
[|(LaO) - (AIO2) - |(LaO)] (or [|(A102) - (LaO) - |(A102)]) as a unit cell, as shown in 
Fig. [TlT c). This is analogous to the unit cell in a type 2 model in Tasker's classification. 
From group theory, it is known that spontaneous polarization can be observed only in 
a direction where the crystal does not have mirror symmetry. Since cubic perovskites 
do have mirror symmetry in the [001] direction (we neglect surface or interface induced 
lattice polarization for now), it is useful to take dipole-free unit cells to reflect the lack 
of polarization in the bulk. 

When we take this unit cell, the stability of a polar surface or interface can be discussed 
by only considering the charge neutrality of each unit cell, since now there is no net 
dipole created by the stacking of charged layers. For example, the instability of the 
(AIO2) terminated surface of LaAlOa is naturally derived because the top-most unit cell 
is [(A102)~ |(LaO)"'"]^°'^, which violates charge neutrality, as shown in Fig. [TU c). Thus 
the (001) surface of LaAlOs cannot keep the bulk termination, either by an AIO2 or LaO 
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49|. Once charge neutrality 



layer, and must reconstruct to compensate this charge 
of all the unit cells is achieved, the system is free to undergo electron / hole modulation 
or interdiffusion / displacement of the atoms, which creates only finite dipoles and does 
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Figure 12: Schematic illustrations of the charge structure across the two types of bulk- 
frozen LaAlOa/SrTiOs interfaces, assuming the ionic charges based on the valence states 
in bulk LaAlOs and SrTiOs. Taking dipole-free unit cells, the interface unit cell cannot 
keep charge neutrality, and the simplest neutral interface stoichiometry is written in the 
right, (a) LaO/Ti02 terminated LaAlOs/SrTiOa interface and (b) A102/SrO terminated 
LaAlOs/SrTiOa interface. 



not violate neutrality^. These perturbations are important because they determine the 
real charge structure, for example via lattice distortion close to the surface of LaAlOa 
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49|. 



Following these arguments, the polar nature of given bulk frozen surfaces and inter- 
faces is clearly defined, considering local charge neutrality using dipole-free unit cells: a 
bulk frozen surface or interface is polar if the unit cell at the surface or interface cannot 
keep charge neutrality when we take dipole-free unit cells in the bulk. 



4.2 LaAlOa/SrTiOs in the local charge neutrality picture 

Taking dipole-free unit cells of perovskites ABO^ in the [001] direction, namely [|(j40) 
- (-BO2) - ^(AO)], the interface unit cell can be treated as a (5-dopant at the interface. 
As shown in Fig. [12], the interface unit cell does not keep the stoichiometry of either 
of the bulk materials, not even a simple mixture of them. This issue is actually one 
of the central opportunities of the interface science of heterostructures. Namely, the 
LaO/Ti02 terminated LaAlOs/SrTiOs interface has Lao.sSro.sTiOa as the interface unit 
cell. Considering the formal electronic charges of La'^"'', Sr^"*", and 0^~, the Ti ion in this 

unit cell should take Ti^'^"*" as a formal valence, which is the same reconstructed state as 

^Thc effect of diffusion is further discussed in Section 16.11 as an example of sources of such finite 
dipoles. Another source of interface dipoles, the quadrupolar discontinuity, is discussed in Section [6.31 
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- Ev ^ Ev 

SrTiOs LaAIOs 
Lao.sSro.sTiOs 

Figure 13: (a) Schematic band structures of SrTiOs, LaAIOs and Lao.sSro.sTiOa, as- 
suming that SrTiOa and LaAlOa are intrinsic semiconductors, (b) Schematic energy 
band diagram of the LaAlOs/SrTiOs interface with the (5-dopant Lao.sSro.sTiOs at the 
interface. 



that for the previous discussion by the polar discontinuity picture. 

Similarly, the A102/SrO terminated interface has Lao.sSro.sAlOa as the interface unit 
cell, and due to the fixed ionic charges of the elements, such a unit cell is not charge neutral 
and thus unstable: {Lal\Srl'\A\^~^Ol~)^'^~ . Instead, allowing a change of the oxygen 
number in the interface unit cell, Lao.5Sro.5AlO2.75 is charge neutral^°, and the decrease 
of the oxygen content at the interface to stabilize a p-type LaAlOs/SrTiOa interface 
is naturally derived. Note in this case, these oxygen vacancies are introduced to keep 
charge neutrality, and thus do not provide any free electrons. Thus, the electrons/oxygen 
vacancies to reconcile the polar instability of the bulk frozen LaAlOs/SrTiOs interfaces 
are supplied by the interface unit cell itself. 



50|, that is 



In semiconductors, (5-doping is usually achieved in a symmetric geometry 
the dopant layer is sandwiched between the same host material. In the LaAlOa/SrTiOs 
case, by contrast, the two sandwiching materials have different band structure, with 
SrTiOa having the narrower bandgap. Therefore, broadening of the charge distribu- 
tion from the 5-dopant occurs only in the SrTiOa side 461]. Figure [T3] shows how the 



LaAlOa/SrTiOa interface can be depicted in a semiconductor energy band diagram. 



^''This composition is not stable in bulk pcrovskitc form, but can be considered as a mixture of the 
bulk compounds, LaAIOs, Sr3Al206, and SrAl204, stabilized epitaxially at the interface. 



22 




(b) 

vw 



Ec 
Ef 
Ev 



Figure 14: Schematic band diagram of an (a) undepleted and (b) depleted 5-doped semi- 
conductor. The undepleted structure contains free carriers as well as ionized impurities. 
The depleted structure contains the same amount of donor and acceptor impurities (Goss- 
mann and Schubert [5( 



4.3 Coupling of polar discontinuities 

When two polar discontinuities are brought in proximity to one another, coupling of the 
charges can occur to minimize the total energy of the system. This is just like the coupling 
of dopant layers in (5-doped semiconductor heterostructures, which can be understood in 
terms of depleted and undepleted structures. Figure [TW a) shows the band diagram of 
an undepleted semiconductor with one layer of 5-dopant. Since the dopant is positively 
ionized, an equal number of free electrons are left, and they screen out the potential 
created by the (5-dopant. As a consequence, the structure is neutral, and has zero electric 
field sufficiently far away from the 5-dopant layer. 

On the other hand, Fig.fl^b) shows the band diagram of a depleted (5-doped structure, 
where the number of donors is equal to that of acceptors, and they are sufficiently close 
in space. As a result, all the free carriers recombine and the structure is depleted. The 
critical parameters to treat electronic coupling of two 5-doping layers are the distance 
between them and the dielectric constant of the medium. When the distance between 
them is smaller than the length scale of band bending of the undepleted structure, they 
couple and the system goes to a depleted state. Note in this case, charge neutrality in the 
neighborhood of one dopant layer is not necessarily maintained. Superlattice calculations 
using density functional theory show that the above threshold is observed by changing 
the thickness of each layer in LaAlOs/SrTiOs superlattices, which can be captured in 



terms of a simple capacitor model 
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4l|. 



In 5-doped semiconductor heterostructures, only considering modulation of free carri- 
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ers is sufficient to describe the coupling of the (5-dopant layers. However, when the system 
allows possibilities of excess charges other than the free carriers, which come from outside 
of the constructed crystal - e.g., anion or cation vacancies or foreign atoms absorbed to 
the surface, the electrostatic potential can be minimized via them. For example, the 
instability of the polar A102-terminated LaAlOs surface can be solved by introducing 
positively charged surface oxygen vacancies 49|. The LaAlOa thickness dependence of 



LaO-Ti02 terminated LaAlOs/SrTiOs can be explained by a simple assumption, where 
we only consider coupling of the free electrons provided by the LaAlOs/SrTiOs interface 
and the surface oxygen vacancies to keep the local charge neutrality of the polar LaAlOa 
surface. Figure [TS] shows a schematic structure of the LaAlOs/SrTiOa heterojunction, 
where two polar discontinuities exist at the LaAlOa/SrTiOs interface and the LaAlOs 
surface. When the LaAlOa ffim is sufficiently thick [Fig. fTM a)] these polar discontinuities 
are decoupled and charge neutrality is preserved locally by introducing oxygen vacancies 
at the surface and taking the Ti valence of 3.5+ at the interface unit cell. Fractionally 
filled Ti valence provides itinerant electrons, and therefore the system is metallic in the 
thick limit. 

If the two polar discontinuities are brought closer [Fig. [TSlfb)]. the external charges 
and the free carriers (surface oxygen vacancies and extra electrons in Ti valence) can 
recombine in an environment with oxygen gas, and the system is depleted and the con- 
ductivity disappears^^. Here, the word "deplete" is used to mean "reduce the amount 
of charge other than charge bound to the crystal" , and the full depletion of the LaAlOa 
surface means extinguishing the positively charged oxygen deficiency at the surface — 
i.e. the surface turns to an unreconstructed bulk frozen state. Manipulatin g th is tran- 



sition appears to roughly capture the essence of writing nanoscale features 



corresponds to the "writing" of surface charge 



52|. 



4^, which 



^^This is similar to the recombination of free electrons and holes, and the actual recombination can 
be written by following Kroger- Vink notation 5l| as: ^0-2 + 2e' + Vq 0. 



24 



(a) 



| AI01.75| 

LaO 
AIO2 
LaO 
AIO2 



LaO 
AIO2 
LaO 



AIO2 



LaO 



lTi^^"02l 



SrO 



SrO 



(b) 


P 


1 AIO2 1 


-0.5 


LaO 1^ 




1 AIO2 1 





LaO 1^ 




1 AIO2 1 





LaO 1^ 




|Ti''"02| 


+0.5 


1 SrO H 




|Ti''"02| 





1 SrO M 





Figure 15: Schematic illustration of the charges (p) of the dipole-free unit cells of atom- 
ically abrupt LaAlOs/SrTiOs interfaces in (a) the thick LaAlOs limit with the polar 
discontinuities locally neutralized, and (b) the thin LaAlOs limit with depleted polar 
discontinuities. 



4.4 Modulation doping by a proximate polar discontinuity 

Even in systems containing only one polar discontinuity, coupling between the polar 
discontinuity and the other layers can occur in analogy to modulation doping by a S- 



doping layer 53|, as shown in Fig. (TBI At the semiconductor interface, lineup of the 
conduction and valence bands should be maintained, as well as a fixed chemical potential, 
which causes the modulation of carriers resulting in band bending. 

Assume an interface between two intrinsic semiconductors A (narrow bandgap) and B 
(wide bandgap). If the (J-dopant in B is sufficiently far from A, the flat band condition at 
the interface is maintained [Fig. [TBT a)]. When A is brought in proximity to the (5-dopant 
in B, in order to keep the conduction band lineup, the conduction band minimum of A 
lies at lower energy than the chemical potential in B [Fig. [TBT b)]. Since the free carriers 
around the 5-dopant in B have higher energy than the conduction band minimum of A, 
they transfer to A and band bending occurs in A [Fig. HSJ^c)]. As a result, the 5-dopant 
in B is depleted, and A is doped close to the interface. 

An experimental example of such charge modulation was observed from a type 3 polar 
surface of LaAlOs (001), which acts as a (5-dopant, to a narrow layer of the Mott insulator 
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Figure 16: Schematic band diagram of the interface between two intrinsic semiconductors 
A and B, with a layer of (5-dopant in B. (a) A is sufficiently far from the (5-dopant in B. 
(b) Model with no charge modulation though the distance between A and the 5-dopant 
is close, resulting in the mismatch of the chemical potential, (c) Charges transfer from 
B to A, so as to match the chemical potential. 
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Figure 17: (a) Schematic band diagram and crystal structure of a LaVOa quantum well 
embedded close to an A102-terminated LaAlOs (001) surface, (b) Illustrations showing 
how reconstruction charge at the LaAlOs surface is transferred to the buried LaVOs 
quantum well. In order to solve the instability caused by the polar nature of AIO2- 
terminated surface, positive charge is required. When the LaAlOs cap is sufficiently 
thick (left), the LaAlOs surface and the LaVOs well layer are separated, and the positive 
compensating charge remains at the surface. For a thinner spacing (right), the LaVO.s 



well layer accommodates this positive charge, which is energetically more favored |45l. |54 



45 



54| . This system is noteworthy in 



LaVOs with smaller bandg cip, clS shown in Fig. [T7] | 
the discussion of polar discontinuity effects, in that the observed hole-doping can neither 
arise from oxygen vacancies nor inter diffusion. 



4.5 Advantages of the local charge neutrality picture 

As discussed, the local charge neutrality picture gives a clear and self-consistent explana- 
tion for the various phenomena at surfaces and interfaces. The stability of given surfaces 
and interfaces can be simply judged by looking at the local composition. When the in- 
terface composition differs from that of the bulk or a simple superposition of them, a 
different electronic and/or atomic state can be expected. Enforcing local charge neu- 
trality at the interface unit cell naturally explains the source of the charges which can 
change the stoichiometry or the electron number of the constituent atoms from those of 
the bulk. For example, in a LaO/Ti02 terminated LaAlOa/SrTiOs bulk frozen interface, 
the interface unit cell is Lao.sSrg.sTiOs and the Ti is 3.5+ to achieve charge neutrality. 
Therefore, the difference between the Ti^"*" in the bulk SrTiOa and the interface Ti^'^+ 
comes from the interface unit cell itself, and not from anywhere else. 

Another advantage of taking the local charge neutrality picture becomes clear when 
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we discuss the coupling of polar discontinuities. The stability of the system can be 
discussed through the distance between the polar discontinuities and the screening length 
of the host material. In other words, the stability is determined by the balance of the 
activation energy of the dopant and the electrostatic energy allowing polarization of the 
media between the dopant layers. Therefore, the total polarizability of the media can 
be considered in the calculation, and is connected to the bulk permittivity in the thick 
limit. Note when the media is thin, the local effective permittivity is non-trivial since the 
local atomic displacements can be different from that in the bulk, and the local dielectric 

n 

approximation breaks down [55[. 



5 Equivalence of the two pictures 

Thus far we have discussed boundary conditions based on two different choices for the 
unit cell. This was implicitly taking a microscopic viewpoint, since the electric field E 
and the total charge density ptot were connected by Gauss' law: EqV ■ E = ptot, where 
is the vacuum permittitivity. Since ptot is used (hence the atomic-scale stepped or 
sawtooth potentials), different choices for the unit cell were irrelevant so long as global 
charge neutrality was considered. Here we confirm the equivalence of the two pictures 
from the macroscopic electrostatic viewpoint. 

5.1 Gauss' law for infinite crystals 

In media, treating ptot is often quite complicated, which can be simplified by using the 
electric displacement D: 

D = £oE + P, (3) 
where P is the polarization. Then Gauss' law is given by 

V ■ D = pfree, (4) 
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where pfree is the free charge - the part of the macroscopic charge density due to excess 
charge not intrinsic to the medium, which are designated as bound charge Pbound- These 
definitions do not depend on whether the charges are locahzed or itinerant, and are just 
introduced for practical convenience to treat the displacement of the bound charges as the 
dielectric response of the media to the electric field. Since the total charge is conserved 
(ptot = Pfree + Abound), by taking the divergence of Eq. ([3]), 

Pbound = -V - P. (5) 

In infinite crystals, only the divergence of E, P, and D has physical meaning, and the 
polarization arising from the density of unit cell dipole moments Pdipoic can be neglected 
since it merely adds a constant value to P. 



5.2 Gauss' law for finite crystals 

When the crystal is finite, Pdipoic drops to zero at the surface. Thus the choice of unit cell 
is important, since it changes the nature of the discontinuity at the surface. According 
to Eq. ([5]) the magnitude of the polarization discontinuity at the surface is determined 
by Pbound there. Since ptot at the surface does not depend on the choice of unit cell, 
uncertainty in the dipole moment of the unit cell is compensated by whether the charges 
at the surface are defined to be free or bound. — the surface charge is bound when it 
belongs to a bulk unit cell, and free when not, and this definition does not depend on 
the origin of the charges 12 1 . 



For simplicity, let us adopt the bulk dielectric constant e to connect E and the induced 
polarization Pjnd by the field, namely 

Pind = {e- £o)E, (6) 

hence 

P = Pind + Pdipoic- (7) 
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Figure 18: Schematic illustration of macroscopic D, Pdipole, P, E, pbound, Pfree, and ptot 
close to the (001) surface of LaAlOs. (a) Taking [(A102)~ (LaO)"*"] as a unit cell and (b) 
taking [|(LaO) - (AIO2) - |(LaO)] as a unit cell. 



Substituting them in Eq. ([3]) we obtain 



D = eE + P 



dipole- 



(8) 



Now, let us calculate the value of E, Pdipoic, D, and the charge density, considering the 
unreconstructed A102-terminated (001) surface of LaAlOa as an example, for the two 
different unit cell assignments previously discussed. 



5.2.1 Taking [(AIO2) (LaO)^] as a unit cell — the polar discontinuity picture 

The z-component of the vectors D, E, and P are denoted as D, E, and P. As shown 
in Fig. [TST a). all the ionic charges are included in the unit cells, and thus bound. To 
fix the constant in E, the vacuum can be taken as the reference for i? = 0. Due to 
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the absence of free charges, D = from Eq. (jlj). In this case, the [(A102)~ (LaO)"*^] 
unit cell has a dipole moment of where go is the elementary charge and a is the 
pseudocubic lattice constant of LaAlOa. Considering the unit cell volume a^, Fdipoie is 
given by Fdipoic = -^^(-2)) where 6(z) is the step function. Equation ([H]) immediately 
provides E as a. function of the position, namely E = — ^ 9{z). The total polarization 

is given by P = ( (e — Eq) ■ — ^ + ) 9(z) = —-^6(z), following equations (E]) and 
\ zea'^ 2a^ J e 2a'^ 

(17]), and from Eq. (jSl), Pbound = P = — — -^^{z)^ where 5{z) is the Dirac 5 function. 

az e 2a^ 

This indicates that the system has bound charges of — at the surface. 

€ 2a^ 



5.2.2 Taking [|(LaO) - (AIO2) - |(LaO)] as a unit cell - the local charge 
neutrality picture 

When taking dipole-free unit cells [Fig. fTSTb)]. the topmost unit cell is [(AIO2) - |(LaO)], 



unit cell is free, and pfrcc = o^i^)- By integrating Eq. (jll) and using the boundary 

2a^ 

condition that D = in the vacuum, we obtain D = —-^^6{z). Since Pdipoic = 0, i? is 



which is not a bulk unit cell. Therefore, the half negative charge which belongs to this 

%_ 

%_ 
D P 

obtained as £" = dipoie _ — ^I^q/^\ According to Eq. ([3]), the total polarization 

is P = D — EqE = - — —-^9(z), and Pbound is obtained by Eq. (15]), namely Pbound = 

e 2a^ 

-^^6{z). This bound charge density appeared as a response to the electric field E 

arising from pfrce- Thus we confirm that ptotai is independent of the choice of unit cell. 



6 Further discussions 
6.1 Effect of interdiffusion 

In real systems, interdiffusion of atoms across the interface is inevitable, and we should 
note its effect on the electrostatic stability of the system. Interdiffusion is a process 
where atoms are locally exchanged, and does not change the charge neutrality conditions 
around the interface, except for the finite dipoles induced by the modulation of charges. 
Since the instability of polar interfaces can be derived from the lack of charge neutrality 
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Figure 19: Schematic charge structure p and electrostatic potential of (a) an abrupt 
interface and (b) an interface with interdiffusion. Five dipole-free unit cells are taken 
as an example of a cluster covering the interface region. The dashed line in (b) is the 
duplication of in (a), (c) Difference p(b) — p{a) of the charge distributions in (a) and 
(b), and created potential shift 0dip by the dipole indicated by the shaded arrow. 



around them, simple interdiffusion in a finite region cannot compensate it. That is, local 
stoichiometric interdiffusion can neither create, nor remove, a potential divergence. Here, 
this point is emphasized by considering a simple model. 

Figure. [TW a) shows the charge structure p and the calculated electrostatic potential 
(p of an abrupt polar interface. By taking dipole-free unit cells, the interface unit cell has 
a half positive charge per 2D unit cell, and the system does not keep charge neutrality 
as it is, and is unstable. Figure. [TW b) shows an example of interdiffusion, where half of 
the negatively charged layer and half of the charge neutral layer close to the interface are 
swapped, compared to the model in Fig. [TOTa). In this case, the extra positive charge 
appears at a different position, but the amount of the charge is conserved, and a similar 
instability in arises as in the abrupt case. 

In order to show that interdiffusion does not change the stability or instability of 
a polar interface, it is useful to consider charge neutrality in a cluster consisting of a 
sufficiently large number of dipole-free unit cells covering the interface region. By taking 
dipole-free unit cells, the electrostatic structure in the bulk can be neglected to determine 
the stability of the system. As an example, let us take a cluster shown in Figs. [TW a) and 
(b). Since the modulation of the charge to achieve the interdiffused model occurs inside 
the cluster, the total amount of the charge in the cluster is conserved, and thus the same 
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instability appears in both cases. 



For comparison of these models, let us focus on the difference of the charge distribu- 
tions p(b) — p(a) in the two models, as shown in Fig. ITOT c). where a dipole with a finite 
size appears. Since (p is calculated by spatially integrating the charge distribution twice, 
the difference of in Figs. ITW a) and (b) should be equal to the potential shift created by 
the dipole in Fig. [TW c). Indeed, (j) in Fig- fTW b) has a shift from that in Fig. [TW a). which 
is the same amount as the dipole shift in Fig. [TW c). It should be highlighted that this 
dipole shift can indeed be an interface-specific additional driving force for interdiffusion. 



and has been suggested to fundamentally limit the abruptness of some interfaces 
The energy associated with a band offset, for example, can be reduced by forming this 
dipole. 



Finally, we should note the difference between simple interdiffusion and change of 
the interface composition. In the interdiffusion process, only exchanging atoms in the 
finite interface region is allowed, and therefore the total number of atoms of each element 
in the region is conserved. On the other hand, interface composition can be changed, 
for example by inserting other materials, segregation of atoms, or creating vacancies. 
For example, the reconstruction model in Fig. |6]^b) to compensate the instability of a 
polar Ge/GaAs interface cannot be achieved by a simple roughening at the interface: the 
numbers of Ge, Ga, and As atoms are different compared to those in the abrupt model 
[Fig.|6i:a)]. 



In summary, the effect of the interdiffusion appears as an extra interface dipole mo- 
ment at the interface from the viewpoint of electrostatics. However, it does not change 
the total number of charges at the vicinity of the interface, and thus cannot screen the 
charge imbalance at polar interfaces. In order to compensate the instability of a polar 
interface, therefore, introduction of a compositional change or other compensating charge 
is required. 
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6.2 Role of correlation effects 

So far we have tried to simply describe the perovskite polar discontinuity using semi- 
conductor language by taking dipole-free unit cells and treating the interface unit cell 
as a (5-dopant. It should be mentioned that in the presence of strong electron-electron 
correlations, commonly found in transition metal oxides, it is formally impossible to draw 
semiconductor energy band diagrams based on the single-particle picture of independent 
electrons^^ js?!. Issl]. Also, in order to draw band diagrams, we should know which part of 
the charges are to be assigned as free carriers, which is non-trivial in correlated systems, 
such as Mott insulators. For example, a perovskite with 1 electron per unit cell can give 
rise to an effective carrier density ranging from ~ 1 x 10^^ cm"'^ to zero, depnding on the 
correlation strength. 

However, the local charge neutrality picture can still provide important information on 
the interface charge structure of transition metal oxides, since correlations cannot change 

^may be significantly 



the amount of total charge. While the distribution of this chare 
modified by correlation features in the electronic compressibility 59| , the charge can still 
be determined in a cluster consisting of a sufficiently large number of dipole-free unit 
cells. This is because the material outside of the charge modulation region consists of 
dipole-free unit cells of bulk, which are charge neutral and thus stable. 

6.3 Quadrupolar discontinuity 

Thus far, we considered the stability of polar discontinuities, and showed that the amount 
of charge needed at the interface can be determined by taking dipole-free unit cells and 
considering local charge neutrality in each unit cell. Of course, local charge modulation 
is allowed once neutrality is obtained, since it does not violate global charge neutrality. 
This dipole energy arising from the modulation of charges determines the real charge 
distribution in the system. 

Here, we note that there is an intrinsic dipole shift at the interface arising from the 



^^Owing to the relatively small (still non-negligible) correlation effects in SrTiOa with the Ti 3dP con- 
figuration and weak 2p-3d hybridization in the coherent state of doped SrTiOs [56[ , the schematic band 
diagram shown in Fig.[T3]is still reasonably valid, approximating SrTiOa to be a band semiconductor. 
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Figure 20: (a) Schematic structure, charge density p, and electrostatic potential of a 
LaAlOs/SrTiOa interface with a Ti^'^"'"02 interface layer, showing a finite shift A. Taking 
(b) [(iLaO)+i/2.A102^-(|LaO)+i/2] and (c) [(iA102)-i/2-LaO+i-(iA102)-^/2] as a unit 
cell. The blue and pink shaded areas in (b) and (c) show the unit cells with opposite 
signs of quadrupole moment, respectively, and the green shaded area in (b) shows the 
unit cell with a finite dipole moment, (d) An example of dipole- and quadrupole-free 
unit cells. 



charge stacking sequence, which should be distinguished from this charge modulation. 
For example, consider a LaAlOs/SrTiOa (001) interface with Ti^'^"*" at the interface to 
solve the instability of the polar interface, as shown in Fig. [20^ a). Although there is no 
potential divergence, a finite shift A remains between the two materials when considering 
the averaged electrostatic potential on both sides. 

The origin of this potential shift can be understood based on the discontinuity of the 
quadrupole moment of the unit cells. In one dimension, the quadrupole moment density 
Q is defined as = -P, where P is the dipole moment density. It has the same form as 
that of Gauss' law connecting the dipole moment and the charge — i.e., the source of the 
quadrupole moment is the dipole moment. Therefore, following the same argument as 
in Section 14. ![ a finite interface dipole moment appears at a quadrupolar discontinuity, 
when the quadrupole moment is different on the two sides. 

Moreover, this quadrupole moment is proportional to the potential shift at that point 
in the absence of free charge, since L> = and thus EqE = —P. Therefore if the unit 
cell has a finite quadrupole moment, that means the averaged potential in the unit cell 
is shifted by the corresponding value. Note the quadrupole moment does not induce a 
potential shift outside of the unit cell, while the dipole moment does change. 

For example, [(iLaO)+i/2.Aio^i-(iLaO)+i/2]^ a dipole-free unit cell of LaAlOs (001) 
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[Fig. [20r b)]. has a finite quadrupole moment, while the charge neutral stacking of SrTiOs 
creates no quadrupole moment in the unit cells^^. This difference of quadrupole moment 
between the two sides adds a finite potential shift to the averaged potential in the unit cell. 
Figures EUT b) and (c) show two ways of taking dipole-free unit cells in cubic perovskites, 
where the sign of the quadrupole moment is opposite. The choice of unit cells cannot 
change the electrostatic potential, and indeed this difference of the quadrupole moment of 
the unit cells is compensated by the dipole moment of the interface unit cell in Fig. [2CT b). 

The ambiguity of the quadrupole moment in the unit cells is reminiscent of the un- 
certainty of the dipole moment, as discussed in Section 14. 1^ and we can take the same 
prescription as we did previously — taking quadrupole-free unit cells. It is more compli- 
cated since these quadrupole-free unit cells should not have a dipole moment at the same 
time, in order to avoid the instability of the polar discontinuity, and Fig. [20r d) shows 
one way to take such unit cells. Note it is always possible to take unit cells which do not 
have either dipole or quadrupole moments in one direction in any bulk crystal^^. Once 
we take this dipole- and quadrupole-free unit cell, the size of the interface dipoles can be 
determined locally by considering the dipole moment in the interface unit cell(s), since 
there is no shift in the averaged potential in the bulk unit cells due to the absence of 
quadrupole moment. 



7 Summary 

As pointed out by Tasker, the order of the ionic charge stacking plays an important role for 
the stability of the surfaces and interfaces, and sometimes they show a diverging potential 
when following bulk compositions and electronic states. In semiconductors where the 
number of ionic charges are fixed, they are usually stabilized by atomic reconstructions, as 
shown by Harrison et ai. However, in transition metal oxides, electronic reconstructions 
provide another possibility to reconcile the instability. Exploiting this degree of freedom 

^■^ Absence of covalcncy is assumed. 

^'*This can be easily proved from the periodicity of the charge structure and the macroscopic charge 
neutrahty of the lattice by using the intermediate-value theorem. 
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has been a central topic of recent research in oxide heterostructures, but with much 
debate over the relative contribution of this effect, as compared to those common to all 
interfaces such as defects and diffusion. 

We hope the reader will agree that the attribution of various experimental measure- 
ments in real materials to purely the polar discontinuity, as opposed to chemical defects, 
is not really a valid separation. It is rather governed by the response of the total system 
to the electro-chemical potential (in the chemistry sense), subject to specific thermody- 
namic, kinetic, and geometric boundary conditions. Given the freedom to define the unit 
cell, taking dipole-free units provides a convenient format for this perspective. While 
not explicitly discussed here, interface screening by local lattice polarization can be di- 
rectly incorporated, including subtleties that arise at very short length scales. Thus these 
topics connect to current research in surface and interface phenomena for ferroelectrics 
and multi-ferroics. Furthermore, recent discussion of the role of film stoichiometry can 



601 1 . Ultimately, the requirement to 



be understood together with these other aspects 
stabilize polar discontinuities is independent of whether the charge is free or bound, itin- 
erant or localized, and electronic or ionic. Perhaps the central message of this simple 
review chapter is that one cannot claim atomically precise oxide heterostuctures without 
understanding how the electrostatic boundary conditions are accomodated. This is both 
a problem to be solved, as well as a fascinating synthetic opportunity. 
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